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^ ■ Abstract 
O 

^ ' We investigate (2+l)-d Hamiltonian lattice gauge theory using a class of 

m ' 

Q^ ' Hamiltonians having exactly known vacuum states. These theories are shown 



to have a wide range of possible classical continuum limits which differ from 
that of the standard Kogut-Susskind Hamiltonian. This conclusion is at vari- 
ance with some previously published results. We examine the quantum con- 



^ ' tinuum behavior of these theories by both analytic and numerical methods 



including plaquette space integration and standard Monte Carlo techniques. 
String tension and variational estimates for the J^^ = 0^^ glueball spectra 
are presented for SU(3). We find that in spite of the wide range of classical 
behavior predicted, these theories correspond to only two distinct quantum 
systems in the weak coupling limit. One of these quantum limits gives string 
tensions and glueball states which show scaling in weak coupling which agrees 
with the perturbative prediction for the (2+l)-d problem. 
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I. INTRODUCTION 



Our current understanding of the strong force gives us a picture of fermion fields in- 
teracting via gauge fields generated by an underlying local SU (3) color symmetry. Since the 
physical hadrons are observed to be color singlets, this picture requires exact color confine- 
ment. We also know that strong interaction processes involving large momentum transfer are 
characterized by nearly free hadronic constituents. The current belief is that non-Abelian 
gauge theories such as SU(3)coior simultaneously possess the features of asymptotic freedom 
for short distance phenomena and exact color confinement. Clarification of the low energy 
behavior of SU(3)coior gauge theory requires a nonperturbative approach and is surely one 
of the most fundamental problems in the field of strong interaction physics. 

Truly accurate lattice calculations of glueball masses without dynamical fermions is an 
essential first step toward the long range goal of describing completely the hadronic spectrum 
from QCD. Wilson's nonperturbative formulation of gauge theory on a space-time lattice 
was developed as a means of controlling ultraviolet divergences inherent in QCD. Truncation 
of the lattice then yields a theory having a finite number of degrees of freedom making direct 
computer simulation of the system feasible. Soon afterward, an alternative Hamiltonian 
formulation of the problem was given by Kogut and Susskind In this canonical form of 
the theory, physical observables such as glueball masses and string tension appear directly 
as eigenvalues, and in addition, the gluon degrees of freedom reside on the links of a three 
dimensional space lattice instead of the four dimensional Euclidean space-time lattice used 
in Wilson's approach. 

A potentially important feature of the lattice Hamiltonian is that it may not be unique. 
The common wisdom is that it must be locally gauge invariant and possess a classical contin- 
uum limit in accord with the continuum gauge theory it approximates. Beyond these rather 
general constraints, the precise form of the Hamiltonian is not determined, and this freedom 
can be used to fix the vacuum wavefunction, resulting in a theory having a Hamiltonian 
modified relative to the standard Kogut- Susskind (KS) form This is the approach 
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we have adopted here. It represents a potentially important advance in the Hamiltonian 
formulation of lattice gauge theory since it avoids the inherent computational complexity 
and limited numerical accuracy of either dynamically evolving the exact vacuum state using 
Monte Carlo techniques 0, or introducing into the formalism an approximate variational 
ansatz for the vacuum |§] [0] |0 • With an exact vacuum in hand, variational estimates 
of the glueball mass spectra become rigorous upper bounds because the vacuum energy is 
precisely known. 

In this paper, we investigate a class of (2+l)-dimensional lattice Hamiltonians for the 
SU(N) gauge theory which have an independent-plaquette vacuum structure. Achieving a 
theory wherein the vacuum is exact involves adding a term, AH , to the KS Hamiltonian. 
The notation here is suggestive of the hope that this modification is small in some sense and 
leads to a theory consistent with the corresponding continuum gauge theory. The fact that 
the exact vacuum has a simple structure with no long range magnetic correlations has been 
the source of concern to some authors |0| however, we believe that equally compelling 
arguments can be made in favor of the present formalism 

We view the classical limit of our lattice Hamiltonian as an indicator that perhaps the 
present theory can agree with continuum QCD. In the appendix, we demonstrate that AH 
remains finite in the classical continuum limit of the (2+l)-d theory. One also finds that 
this term has a vanishing limit in the more physically relevant (3+l)-d version of the theory. 
This is encouraging since, after all, we think of the KS Hamiltonian as a valid discrete 
representation of the continuum theory based largely on this same classical correspondence. 
We keep in mind, however, that we are ultimately interested in clarifying the quantum 
behavior of these systems particularly in the scaling region. 

Also in support of the present approach, consider the key role played by the independent- 
plaquette component of the vacuum in the KS theory. Several variational calculations have 
been reported which clearly support the notion of a vacuum dominated by long range mag- 
netic disorder even well into the scaling region 10 [0] 0- While this simple picture 
can not be correct in the extreme weak coupling limit, it does seem to capture the main 

3 



features of the true KS vacuum in regions where scahng is observed and is thought to be 



everywhere exact as a long wavelength representation |TT|] ||T2| |T3|. Again, this allows one 
to argue that at least in the near-weak coupling regime, AH represents a small adjustment 
of the KS dynamics which ensures an independent-plaquette vacuum. 

In this paper, we present numerical results for the SU(3) version of the theory using 
both exact (expectation values of gauge invariant operators are exactly integrable in the 
(2+l)-d theory due to the simple Jacobian) and standard Monte Carlo methods. We make 
use of Monte Carlo mehtods since we are ultimately interested in the (3+l)-d version of the 
SU(3) theory (where a complicated Jacobian involving Bianchi identities appears) and wish 
to calibrate our computer codes where exact results are possible. 

To further motivate the present approach, we have also obtained the SU(2) string tension 
using our formalism which we will now compare with previously published results obtained 



by more conventional means. Teper [|14[ [|T^] has recently obtained = 1.336 ± 0.01 (see 



reference for notation) for SU(2) string tension within the D = 3 Euclidean theory. This 
compares very well with our result of 1.30 ±0.1 (approximate error estimate) for this same 
quantity. Using a one and two-plaquette varational ansatz for the KS vacuum, Arisue et al 
obtained j3y/a = 1.61 for the independent-plaquette vacuum and 1.38 for the two-plaquette 
vacuum. These values are eyeball estimates from their figure 5 in the scaling region near 
(3 = 1.86. Comparision of our result with these variational results clearly suggests that our 
exact local vacuum formalism is doing more than a one-plaquette ansatz does within the KS 
theory. We are very encouraged by this result and proceed to further assess the strengths 
and weaknesses of this potentially important approach. 

The paper is organized as follows. In section P and the appendix, we discuss the classical 
continuum limits of a family of lattice Hamiltonians having exactly known vacuum states. 
Section ^ discusses string tension for the present theory by examining exact results for 
the plaquette expectation value. In section 0, we describe the variational approach to the 
0"'"+ glueball spectra and in section |V| we give the specifics of the Monte Carlo simulations. 



An examination of our numerical results appears in section ^ and finally, section |VI1 



summarizes the main points and discuss prospects for future investigation. 



II. LATTICE HAMILTONIANS WITH EXACT VACUUM STATES 

The Kogut-Susskind form of the lattice Hamihonian for SU(N) gauge theory can be 
derived in a straightforward way from Wilson's lattice action using the standard canonical 



prescription for obtaining a Hamiltonian from a Lagrangian |T^. The result is, 



i/KS = (-E^Et + E[2A^ - Tr(^p + Ul)] . (1) 
zu ay p 

Here, g is the lattice coupling constant, a is the lattice spacing, / labels the links upon which 
gauge degrees of freedom (SU(N) rotations) reside and serve to join neighboring lattice sites, 
a is a color label, Ef are color-electric field operators representing the variables conguate to 
the gauge (link) variables, p labels the plaquettes, and finally. Up are plaquette variables. A 
plaquette is defined as the path-ordered product of link variables, Ui, obtained by traversing 
a closed square loop in the lattice having one link to a side. Sums over repeated indices are 
assumed throughout this paper. 

Being conjugate variables, the link degrees of freedom and the color-electric field opera- 
tors satisfy equal time commutation relations with respect to each other, i.e., 

[Er, U,] = Su'^Ui and [E^, Uj,] = Sn^Ul^ . (2) 

We can see that E" acts as a differential operator on the variable Ui with respect to the 
parameters of the SU(N) rotation. It follows then that the first term of Eq. (|1]) is a sum of 
the kinetic energies possessed by the link variables, and the second is a sum over potential 
terms which depend on the relative orientation of the gauge degrees of freedom through 
couplings provided by the plaquettes. 

In the limit of small lattice spacing, i^KS agrees with the classical Hamiltonian for the 
continuum SU(N) gauge theory. In 3 space dimensions the result is very familiar, 

Hks^^I • + 5° ■ B'')d^x . (3) 



Agreement of the classical limit is usually regarded as a reasonable minimum condition that 
a lattice approximation should satisfy. However, we are ultimately interested in whether the 
quantum dynamics generated by various lattice approximations such as Eq. ([1|) can yield 
the known results from the continuum theory. Thus, we look to the scaling behavior of 
the string tension and glueball masses emerging in the weak coupling limit of our lattice 
calculations for a more meaningful indication of whether a particular lattice formulation is 
capable of producing physical results or not. 

To what extent then is the classical continuum limit indicative of the quantum continuum 
behavior? This is the question of universality in lattice gauge theory. We will demonstrate 
that a particular class of lattice Hamiltonian whose members have a wide variety of possible 
classical limits, correspond to quantum systems which have only two distinct weak coupling 
limits. One of these weak coupling limits appears to be the physical one and is in agreement 
with perturbative predictions from the continuum gauge theory. 

The Hamiltonians we will examine have the desirable feature that their vacuum states 
are exact for all values of the lattice coupling, thereby making the calculation of vacuum 
expectation values of various lattice operators, which are needed for the evaluation of string 
tension and variational estimates of glueball masses, a relatively easy task. An exact vacuum 
can be achieved by adding a term. A//, to the basic KS Hamiltonian, 

H = Hks + ^H . (4) 

At this point, it is convenient to redefine the energy scale by subtracting the constant, 
ultraviolet divergent part of the magnetic term in Eq. (|ID. Also, in (2+1) dimensions, one 
finds that the lattice coupling is related to the dimensionless gauge coupling, e, by = ae^. 
We then have, 

H = '^Ersr --^,1: MU, + f/t) + AH . (5) 

We will choose AH such that the modified Hamiltonian is gauge invariant, Hermitian, and 
positive definite. 
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We will now look for a form for Ai/ which renders a particular vacuum state exact. 
Consider a vacuum of the form, 

l^o) = e«|0) , (6) 

where |0) is the product of individual link ground states defined by E"\0) = 0. In this paper, 
the form we consider for R is, 

R = J2{(3,A, + f32Al), (7) 
p 

with 

A, = Tr([/p + UD . (8) 

Here, (3i and [^2 are as yet undetermined functions of the lattice coupling. If H is positive 
definite and -f^|\l/o) = 0|\l/o), then |\l/o) is assured of being the vacuum. Using this condition 
to determine AH, one finds, 

Ai:^|*o) = [-'^(E-ErR) - ^(ErRmR) + -^E^JI^o) . (9) 
Thus, H has the same classical continuum limit as the KS theory if, 

=2 „2 



AH = -jiE^E^R) - '-{ErR){ErR) + ^E^p ^ . (10) 



p 



We view this condition as a guide in searching for theories which yield the same physics as 
the standard KS form, keeping in mind that we are really interested in the quantum behavior 
of the resulting systems. In fact, one of the results of the present work is that theories need 
not share classical limits to give identical physical results. In the appendix to this paper, 
we look for Pi and (32 such that Eq. ( [T0| ) is true. We find that it can not be satisfied for the 
(2+l)-d problem but that AH can at least be made finite by imposing the constraint, 

along with the condition that (32 diverges no faster than a~^. 



132 < 0{a-^) . (12) 

These conditions eliminate divergences in Ai/ as the lattice spacing vanishes. However, as 
we discuss in that section, they are not enough to fix the classical limit of Aif . One must 
further specify how /5i and (32 individually behave as the lattice coupling and hence the 
lattice spacing are taken to zero. If we relax the latter condition and only enforce Eq. ([TT]), 
then we have shown in the appendix that the resulting Hamiltonians may differ from the 
KS theory by a divergent amount. 

In the next section, we will numerically examine some theories which satisfy Eqs. ( |TlD 



1 /■ 1-1- 1 r 

^ - / • + —- 5" ■ B'')(fx + — ^ / (Px[{y X S") ■ (V X B")] . (13) 



and (|T^ . Specifically, we put /5i = 2c^^ + b and [32 = 4^ [ 20^34 ~ ^] ^'^^ consider the 
continuum limits for fixed values of X and h. Figure |I| indicates the relationships among 
the variables we have introduced thus far for the case 6 = 0. In the appendix, we find the 
classical limit of our modified Hamiltonian to be independent of 6, 
1 1 + ?>X ,0 1 

So, the magnetic component of the resulting Hamiltonian is modified relative to the KS form 
by an amount which depends on the relative strengths of the [3i and (32 terms through the 
parameter X . In particular, the critical value Xc = —\ marks a sign change in the magnetic 
contribution to the Hamiltonian. We also see a contribution which is independent of the /3's 
so long as Eq. (pH]) is satisfied. 

In spite of the widely diverse classical continuum behaviors for these theories, we will 
now argue (and later demonstrate numerically) that Eq. (|Tl|) alone is sufficient to ensure 
only two possible weak coupling limits for the corresponding quantum lattice systems. 

Consider the continuum limit of the vacuum state density \E'o = e^^. This limit will 
determine the continuum behavior of observables such as the string tension as the lattice 
cutoff is removed. From Eqs. (|^) and (0), we can see that as the lattice coupling is taken 
to zero, the extrema of R will dominate integrals over the plaquette degrees of freedom 
(stationary phase condition). In the case of SU(3), there are three turning points for R 
corresponding to Ap = Ap = —2, and Ap = —3. Which one 'wins' in the limit of 



weak coupling simply depends on which density is greater, "^li^p = 6), "^KAp = —2), or 
"^KAp = —3). One can easily verify using Eq. (|16D below that the extreme values, Ap = 6 
or Ap = —3, always dominate the integration. Thus there will be only one phase change 
occuring at the crossover point for these two solutions, 

R(Ap = 6) = R{Ap = -3) 

A + iV/52 = , (14) 

or in terms of X and b, 

A: + 2bCNg' = ~, (15) 

where the critical value Xc = is given by the limit of weak coupling {g 0) and 
corresponds to the sign change in the magnetic part of the Hamiltonian seen in Eq. (|13|). 
We see then that separates two distinct regions of quantum behavior. For example, when 
X < the vacuum expectation value {Ap)o approaches —3 as the lattice approximation 
is removed whereas X > X^. gives a limit in agreement with the classical value of +6 (see 
Eq. ( [A6D ) . In the next section, we present numerical results which are in complete agreement 
with the above analysis. 

So, here is the situation. Each value of X gives a family of lattice Hamiltonians, cor- 
responding to different values of 6, which share a unique classical limit as well as a unique 
limit for the associated quantum system. This supports the concept of universality in the 
sense that all Hamiltonians we have studied which share the same classical limit also yield 
identical physics. However, the mapping between classical and quantum limits is not one- 
to-one. Classical theories corresponding to different values of X map onto only two distinct 
quantum systems in the weak coupling limit. The quantum limit obtained in a particular 
case corresponds directly to the sign of the magnetic contribution to the Hamiltonian in 
the classical limit. Thus we see that many different classical Hamiltonians correspond to 
the same quantum system. The relative strength of the electric and magnetic parts of the 
classical Hamiltonian plays no part in the continuum limit of the quantum system; only the 
sign matters. 
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We now express Eq. (0) in terms of the parameters X and b to be used in the subsequent 
sections. Also, for the remainder of this paper, we will work exclusively with SU(3) gauge 
theory. If we define the lattice coupling parameter, ^ = 

R = ^[^('^ + ^"^p) + - 1^^^)]^^ • (16) 

Before proceeding with the numerical results, we should point out differences between 
conclusions reached in this section and the appendix to this paper and those arrived at 
previously by other authors working with the same Hamiltonians and proposed vacuum 
states used here. The condition given by Eq. (|ll]) is not the same as either of the two 
conditions suggested in the work of Ref. except for the special case of /?2=0. In this special 
case their conditions on a2 (our given by their Eqs. 2.16 and 2.17 are inconsistent, i.e., 
tt2 = lon^ 4^ a2 = respectively. We believe their constraints to be erroneous except 
for their special case of Ri = ai J^p^p with ai = j2c^^- Indeed in this latter case we do 
reproduce their calculations for the mass gap (their Fig. 8) and string tension (their Fig. 
9). The correct treatment of the classical continuum limit is presented in the appendix and 
is well verified numerically here in the results presented in the next section. 



III. SU(3) STRING TENSION 

We are interested in evaluating the quantity. 

Since the Jacobian J' for the change of variables {dUi} J{dUp\ is unity for the (2+l)-d 
problem and the integrand is explicitly a function of Up, the expectation value above reduces 
to an integral over a single plaquette variable, 

J P 

where we have used Eq. (|16]). For SU(3), these two integrals can be reduced to two dimen- 
sions using Up in diagonal form and evaluated directly for various fixed values of X, b, and 

10 



Figure ^ shows a portion of the surface {Ap)o{X,^) with 6 = and clearly shows an 
approach to a discontinuity across = — | in the weak coupling limit as suggested in the 
discussion leading to Eq. (p!5|). For X > X^, the theory is consistent with the classically 
predicted limit of +6, while for X < Xc, the limit —3 is seen. Similar results occur for 
non-zero values of b as pictured in Figs. |^ and ^ where the results for b = ±1 are pictured. 
Again, A'^ = — | emerges as a critical point separating two domains of convergence in the 
limit of weak coupling. 

In (2+1) dimensions, the string tension as obtained from the area law behavior of Wilson 
loops can be shown to be a function of {Ap)o alone [T^. A sufficient condition for this 



result to hold is that the vacuum factorizes into pieces associated with each plaquette as it 
clearly does in this case. One finds that the expectation value for a loop enclosing an area 
m measured in units of is given by the m^^ power of the fundamental plaquette. We then 
have, 

a^K = -^ln(^) . (19) 

The factor | appears because the area law coming from spacelike Wilson loops corresponds to 
an octet (adjoint representation) color source and sink rather than the triplet (fundamental 
representation) q - q potential that we really want ||18|| . 

In addition to the approach described above, we employ a more general method of 
extracting the area law behavior which is also applicable to the more complex (3+l)-d 
problem. Here we calculate the string tension directly from the expectation values of a 
sequence of increasingly larger square Wilson loops. When the loops are large (compared 
to a correlation length on the lattice), we expect the behavior to be dominated by area and 
perimeter laws, 

0/ ^ exp(-A - ^Ka^P - Ca2I) , (20) 

where / indicates an JxJ loop and for fixed values of the lattice coupling. A, K, and C are 
constants. The factor | appears as in Eq. ([191) . So, a quadratic least squares fit to the data, 

11 



- ln{(j)i) + ^Ka^I^ + Ca2I , (21) 

gives A, \d^K, and 2aC for various values of ^. 

The asymptotic scahng behavior obtained from a perturbative expansion of the renor- 



malization group function for the SU(3) theory in 2+1 dimensions is = constant [T^. So, 
we will look for constant weak coupling behavior for lattice quantities such as string tension 
and glueball masses which are measured in units of the lattice spacing as an indication of 
consistency with the continuum theory. 

Thus, we can see just from the results in Figs. and Eq. ([T9|) that only the theories 
corresponding to X > Xc give meaningful string tension results. 

We evaluate the integrals for the string tension calculation suggested by Eq. (^) using a 
Monte Carlo simulation on a finite 2-d space lattice. In the results section, we then compare 
string tension obtained in this way with the 'exact' results as a test of our implementation 
of the statistical integration technique and of the quadratic fitting approach to the string 
tension. This is of interest since in a subsequent paper, we will report results for the (3+l)-d 
problem where comparisons with exact results are not possible as they are here. 



IV. VARIATIONAL ESTIMATE OF SCALAR GLUEBALL SPECTRUM 

Since the eigenvalues of H in the continuum limit correspond to physical observables 
only when gauge invariant states are used, we begin by projecting a set of gauge invariant 
basis states orthogonal to the vacuum state, 

\^^) = (0, - (0i)o)|*o) , where {(Pi)o = (^o|0i|^o) ■ (22) 

It is easy to demonstrate that the trace of any product of link variables constructed by 
traversing a closed path in the lattice is gauge invariant and is thus a potential basis function. 
As there are infinitely many candidates for (pi, an attempt must be made to choose a small 
subset of these which sufficiently spans the space occupied by the low-lying excitations of 
the theory. Obviously, the most efficient choice would be an orthogonal set of basis states. 
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As there appears to be no easy way to choose such a basis beforehand, we will work with 
a nonorthogonal basis and variationally minimize the excited state energies (M) by solving 
the generalized eigenvalue problem, 

det|H-^O|=0, (23) 

r 

with the definition of the Hermitian matrices 

2a 

O,, ^ (^,1^,.) = - . (24) 

Apart from the factor p^, the first matrix is just the Hamiltonian evaluated in the chosen 
basis and the second matrix gives the overlaps among the basis states and is a measure of 
the nonorthogonality of the basis. Thus, for a basis containing n states, there are 'n? + n 
integrals of the form, 

to evaluate. Evaluation of these matrix elements can be simplified in the same way as 
described for {Ap)o in the previous section. However, we use the less efficient Monte Carlo 
approach since we are also interested in the (3+l)-d problem where the integrals can only 
be evaluated statistically. 

The variational basis we choose consists of all square loops having up to 7 hnks to a side. 
They can be expressed as, 

1 

cl>i = jj^Y.^Rei:T{Uu), (26) 

where the lattice has L links to a side and Uu denotes the path-ordered product of link 
variables obtained by traversing the perimeter of an JxJ square loop located at lattice 
position i. Figure ^ gives a graphical representation of these 7 states. In the weak coupling 
(continuum) limit, where the rotational symmetries of continuous space are restored, this 
basis produces estimates of the J^*" = O"*""*" glueball spectrum. 
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V. MONTE CARLO SIMULATION 



Mont Carlo calculations were performed using a 21^ lattice having periodic boundary 
conditions fl^ This lattice size is large enough to give finite size effects for the largest 



loops in our basis which are negligible compared to the statistical fluctuations in our Monte 
Carlo estimates. We determined this by repeating most of the computations presented in 
this paper with smaller 14^ and 7^ lattices. Differences between the 21^ and 14^ results were 
well within statistical noise. It was only with the 7^ lattice that large flnite size effects were 
seen. Thus, we believe that the lattice used here is a very good approximation to an infinite 
one for the basis under consideration. 

We looked at theories corresponding to the four parameter values, X = 0.5, 1.0, 1.5, and 
2.0 with 6 = in each case. For each X we calculated the matrix elements of Eqs. ( |2^ 
for 16 values of the coupling constant distributed over the range 0.30 < ^ < 2.00. We 
employed a multi-hit Metropolis algorithm to generate an ensemble of quasi-independent 
lattice configurations, {Uu} distributed according to the probability density, 

^o(^^) (27) 

During the Metropolis sweeps of the lattice, each link was updated (hit) three times. The 
probability for accepting a link update was maintained near the value 0.5 by dynamically 
adjusting the step (SU(3) rotation) size. Thus, after each sweep, the probability for a link 
to remain unchanged is only 0.5^ = 0.125, increasing the independence of subsequent lattice 
configurations compared to a single hit approach. For each value of ^, thermalization of a 
starting lattice configuration was followed by the generation of an ensemble of 84,000 lattice 
configurations partitioned into 20 subensembles of equal size. Typically, the starting lattice 
was taken from a prior simulation at a nearby value of the coupling in order to reduce the 
number of sweeps required to sufficiently approach the desired distribution. 

Monte Carlo estimates of the expectation values are then simple ensemble averages, 

1 ^ n 

{om)o = -j:om±^, (28) 

i=l V J-^ 
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where = 84000 in our case, and the estimate of the statistical error, ctq, is calculated as 
the variance of estimates obtained from the 20 subensembles. 

VI. RESULTS 

Vacuum expectation values of the square loops in our variational basis as functions of 
the lattice coupling ^ are displayed in Fig. || for A" = 0.5, 1.0, 1.5, and 2.0 with b = 0. 
In the strong coupling limit, where ^ ^ 0, the link variables can be expected to uncouple 
and distribute themselves uniformly over the SU(3) manifold. Our basis states, which are 
products of the links, must obey this same distribution and therefore have zero trace on 
average. The figure shows that the expectation values do approach zero in this limit. In 
weak coupling, where ^ ^ oo and the correlation length diverges (a 0), neighboring 
links should become increasingly correlated. In this limit, the normalized traces of Eq. ( |26D 
approach unity for all /. The figure shows how, as the correlation length increases, so do 
the loop expectation values. Our basis states are becoming more and more alike and the 
off-diagonal elements of the overlap matrix continue to grow in size as the lattice spacing 
shrinks. If ^ were taken large enough, numerical singularity of the overlap matrix would 
eventually occur. 

At a given value of the lattice coupling. Fig. |^ demonstrates that the links become more 
highly correlated as X increases. So the approach to continuum behavior of the system 
(with increasing ^) appears to be more rapid for larger values of X. 

The area law behavior of these 7 states was extracted from the Monte Carlo data both 
by least squares fitting as suggested in Eq. (pTD and from the plaquette expectation value 
as in Eq. (plQl). The resulting string tensions are displayed in Fig. ^ along with the exact 
result. Agreement between Monte Carlo and the exact result is nearly perfect except for 
some significant statistical errors in the least squares result in the strong coupling region. 
Here, the larger loops have large fractional uncertainties (because their numerical values 
are close to zero) and the least squares method requires at least three well determined data 
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points. Some least squares estimates at ^ = 0.30 and 0.40 are missing because the data was 
not accurate enough for the fitting routines to function properly. 

For all four values of X, the resulting string tensions agree both in the strong and weak 
coupling limits. It is only in the transition between the two regimes that the theories differ. 
They appear to be approaching a common constant value of a\fK^ = 0.256 in the continuum 
limit. Of the four X values considered, X = 1.5 shows the strongest appraoch to continuum 
behavior becoming constant near ^ = 0.4. 

Next, we diagonalize our Hamiltonians within the space defined by the variational basis. 
Masses associated with the four lowest eigenstates along with the mass-to-string tension 



ratio are displayed in Figs. § through |IT] for each of the four values of X we considered. 
Errorbars are included in these figures but are nearly invisible since they are small compared 
to the size of the symbols and the vertical scale of the plots. 

The mass gap (lowest eigenstate) enjoys a rather large scaling window extending ap- 
proximately over 0.40 < < 1.4 for X = 2.0 . This is to be contrasted with the result 
for X = 0.5 where scaling set in only for ^ > 0.8 . Again we see that larger X results in 
continuum behavior at smaller values of ^. For ^ > 1.4, basis truncation causes deviation 
from scaling since the decreasing lattice spacing implies that progressively larger loops must 
be included in the variational basis in order to accurately represent the glueballs. The mass- 
to-string tension ratio curves are similar to the mass curves except that the differences in 
the scaling region among the four X values is enhanced by the slower convergence of the 
string tension for small X. 

Results for the second lowest eigenstate are very similar to those obtained for the mass 
gap except that the scaling windows extend only to about ^ = 1.0. This is to be expected 
since excited states should be physically larger than the ground state (as in other familiar 
bound state problems such as the hydrogen atom) and should then suffer from the limitations 
of the basis truncation at smaller values of ^ than in the ground state case. In fact, for 
X = 0.5, the second eigenstate fails to scale at all. An approach to scaling can be seen 
for the third eigenstate in the region 0.6 < ^ < 0.8 for the three X values larger than 0.5. 
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Finally, the fourth eigenstate seems to be beyond what our limited basis can accurately 
represent. 

Examining Figs. ^ and [11] reveals an interesting trade off with regard to choosing the 
'best' value for the parameter X. It is due to the fact that the variational basis is different 
for different values of X. Recall that the set of projectors, 0,, are the same in every case but 
that the basis consists of the parts which are orthogonal to the vacuum which is a function 
of X. So, the question then becomes which X will give the best basis for representing the 
glueballs. Our results show that although larger X produces better scaling in the sense that 
it occurs over a larger window in ^, one appears to pay for this by getting somewhat poorer 
variational estimates of the excited state energies. For example, the fourth eigenstate shows 
approximately a 10% difference between the energies coming from the theories corresponding 
to X = 0.5 and X = 2.0. The lowest eigenstate seems to be unaffected by the modification 
of the basis for different X; that is, the theories agree on the size of the mass gap. 



VII. SUMMARY 

The results obtained here show that there is a class of lattice Hamiltonians for SU(3) 
gauge theory which, in spite of having a range of different classical limits, have the same 
scaling limit for the J^*" = O^"*" glueball masses and string tension. This universality holds 
provided that the vacuum expectation value of a plaquette also has the correct limit, which in 
the present calculations occurs for X > — |. Calculations using plaquette space integration 
2l| show that similar conclusions hold for SU(2) as well as for SU(3) so we anticipate it is 



true for all N. 

Due to an error, previous results in this field did not show true universality in (2+1)- 
dimensions, nor was it realized that there is a continuous class of Hamiltonians (i.e. any value 
of A* > — |) which yields universality. The present work allows us to set the stage for (3+1) 
dimensional lattice calculations with the knowledge that this class of lattice Hamiltonians in 
(3+1) will have the same classical continuum limit as the Kogut-Susskind Hamiltonian. The 
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latter is already known to have the desired classical continuum limit as QCD. Calculations 
are in progress using Monte Carlo techniques in (3+1) dimensions for SU(3) and will be 
reported on later. 



APPENDIX: CLASSICAL CONTINUUM LIMIT 



In this appendix we examine the (2+l)-d classical continuum limit of the vacuum fixing 
term, /\H given in Eq. (p!0|), 

= -iriEtE^R) - i-{ETR){ErR) + -^E^p > (Al) 

where the form of R is defined by Eq. (|^). Using Eqs. (|^), (P), and (||) one can easily verify 
that, 



E^R = Y,{P, + 2(3,A,){ErA,) 

pDl 



(A2) 



and 



E^E-R = E[(/?i + 2P2A,){E^E^A,) + 2P,{E^A,){E^A,)] , (A3) 

where ZlpDZ is a sum restricted to plaquettes sharing the link /. AH can be separated into 
the following three terms, 

A^^i = y E E CnWi + 2P,A,)A, - E 
AH,= ''^Y.i:mE?Ap){EtAp) 

= ^ EE(A + 2P2Ap){E^Ap)r , (A4) 

^ I pDl 

where Cn = ^A^A" = ^ in the expression for AHi is the Casimir invariant. We will now 
take each of these terms in turn and examine the behavior as a — » 0. 
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1. AHi 



Since Trf/p is invariant under cyclic permutations of the four link variables from which 
Up is constructed, the double sum over any function of this trace, J2iJ2p3i /(Trf/p), is just 
four times the sum over plaquettes, 4X]p /(Trf/p). AHi of Eqs. (^) can then be expressed 
as a simple plaquette sum. 



AH, = 2e'CN E(A + 2/32Ap - _^ \ J Ap . (A5) 

p ZL/j^d 6 



The function Ap has a well known continuum limit which is obtained by making a Taylor 
series expansion of the link functions and collecting like powers of the lattice spacing. The 
result is, 

Ap^2N - ^e2a^G?2(^)G?2(^) + 0{a^) . (A6) 

The eight gauge field tensors associated with a plaquette in the 12 plane at lattice location 
X are defined as, 

G^2(f) = diA^ix) - d2A1{x) - er^~<A^,{x)Al{x) . (A7) 

Here, Af{x) denote the two spatial components of the A^^ — 1 gauge potentials, and f'^^"' are 
structure constants defined via commutation relations among the SU(N) group generators, 

[^\]=^r''^■ (A8) 
Using Eq. ( [A^ ) together with X]p o? ^ (Px, the leading terms in AHi are, 

AHi = J (Px{2N - \e'^a^G^^{x)G^^{x) + 0{a^)) 

{Pi + ANI32 - - /^se^a^G- (f)G'?2(x) + Oia^)) . (A9) 

All terms greater than 0{a^) can be eliminated by choosing the parameters /3i and /32 such 
that, 
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which gives, 




(All) 



Notice first the above constraint does not determine the a dependence of /?i and /32. The 



should remain always positive so that the lattice coupling is real. We see then that the 
classical limit of our Hamiltonian is not fixed by Eq. (|A10|) alone; one must further specify 
how f3i and (32 individually depend on a. The possibilities fall into three catagories. First 
if /32 < 0(a~^), then AHi vanishes. This can be realized by choosing P2 =constant, for 
example. Secondly, in the event that (32 > 0(a~^), AHi is divergent for small a. Finally, 
when /?2 is exactly of order a~^, AHi has a finite non-zero limit which depends on the precise 
form chosen. As an example we may put Pi = 2CNa?e'^ ^"^^ -^2 = sAfcwa^e^ which case the 
classical limit for AHi becomes. 



We write the color-magnetic field as a vector to make clear the connection with the more 
familiar case of 3 space dimensions. In (2+l)-d, 3°" has only one component. We see 
that X ^ 1 implies a modification of the magnetic part of the Hamiltonian relative to the 
standard Kogut-Susskind form. Later in this section, we discuss the implications of this 
result. 



Since Tr(£'f Ap) 7^ Ti{EfiAp) in general for I ^ I' G p, we can not directly make the 
replacement J2iJ2p^i ^J2p as we did in the analysis of AHi. Instead, we let Ui,U2,U3, 
and U4 be the four links which define a particular plaquette and define. 



only 'additional' restriction on the form of the /?'s is that the left hand side of Eq. ( |A10[ ) 




(A12) 



2. AH2 



W = UiU2uM 
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Up2 = U2UMU1 

[/p4 = U.UsUM . (A13) 



Now the double sum in AH2 from Eqs. ( |A4| ) can be expressed as an unrestricted sum over 
plaquettes, 

P i=l 

= e^E/32E{Tr[^(^,, - Ul,)]r , (A14) 
p i=i ^ 

where we have also used the definition of Ap and Eqs. (@). When a is small the trace 
appearing in Eq. ( [A14| ) can be expressed, 



iWG^^Tri-— ]+0(a^) 
iiea^G^a + 0{a^) , (A15) 



where the sign ± depends on the sense of rotation in the definition of Upi. In Eqs. ( |A13| ), f/j 



and Up2 involve G12 whereas Up^ and [/p4 circulate in the opposite sense and involve G21 = 
—G12. In the final step above we used the normalization of the generators, Tr(^^) = 
The leading term of Eq. ( |A14|) can now be written, 

AH2 = e^(32j2T.{±^e^'G1,} 

p 1=1 



2 



4e^a'f32 / d'xG'^^{x)G'^2{^) , (A16) 



which comparing to Eq. ( [All| ) is proportional to the limit for AH 



AH2 = --^AHi , as a ^ 

= -]aH^ , SU(3) . (A17) 
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3. AHs 



Finally, consider AH^. In order to decide its limit, let p' and p" be the two plaquettes 
which share a common link /. Furthermore, we define Up' and Up" such that / is the first 
link occuring in the product of four. That is. Up' = UiUmUlUl, and similarly for Up". The 
expression for A//3 appearing in Eqs. ( |^ ) then becomes, 

AHs = ^ ^{(A + 2P2Ap'){ErAp') + {P, + 2P2Ap"){ErAp")y 
^ I 

= J E{(A + '2P2Ap')TT[^{Up' - Ul,)] + (A + 2(32Ap")TT[^{Up" - Ul„)]}' . (A18) 
Using the result of Eq. ( [A15| ) we have, 

Tr[Y(f/p' - Uj,,)] = tea^C^^ix') + ©(a^) . (A19) 

and 

Tr[y - Ul„)] = zea^G^.ix") + ©(a^) . (A20) 

The relative ordering of the tensor indices in the two expressions above (12 vs. 21) occurs 
because p' and p" are defined such that they have opposite rotational sense. AH3 can now 
be written as, 

AH, = + ANP,)[GUx') + G%x")] + ©(a^)}^ 

2 I 

(A21) 

where in the second step we have used Eq. ( |A1CI| ). In the continuum limit, x" = x' + dx 
and the two field tensors appearing in the brackets must have a relative minus sign between 
them from antisymmetry. It follows then that this sum is just the spatial derivative with 
the appropriate power of a. Futhermore, the sum over lattice links can be replaced by 
an area integral together with a sum over the 2 spatial coordinates. That is, we can use 
^ I J2i=x,y The result is, 
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As in the case of AHi, we have written the vector form for B°' even though in the current 
formulation it has only one component. The most significant feature of this result is that 
for any choice of the /5's subject to the constraint given in Eq. ( |A10| ), AH3 makes the same 
contribution to the Hamiltonian. 



4. result 

In (2+1) dimensions, the modified Hamiltonian under consideration has a classical con- 
tinuum limit which differs from the standard Kogut-Susskind limit in the following way. 
Recall that the modified Hamiltonian is given hj H = i^fxs ~ AHi — AH2 — AH^. Based 
on the preceeding analysis of these three terms, the classical limit for this Hamiltonian is 
dependent on the exact form chosen for the /5's. If we decide that we want a finite limit, then 
Eq. ( lAlOl) along with the condition /?2 ^ 0(a~^) must be satisfied. An interesting choice 
which obeys these conditions was given in terms of the parameter X. With this choice we 
find that the magnetic term is modified relative to the Kogut-Susskind form in a way which 
depends on X, 



H '^-J{E''-E'' + — ■ B'^)d'x + Q , (A23) 



where Q is shorthand for the A'-independent result for A//3 given by Eq. ([A22| ). We see that 



the critical value = —\ marks the boundary between two regions where the magnetic 
contributions come in with opposite signs. 
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FIGURES 

FIG. 1. Relationships among variables used to specify the Hamiltonian. To simplify the picture 

we consider here the case 6 = 0. The two darker grey areas indicate regions with distinct weak 
coupling — > oo) limits for the quantum system. For example, X > ^ and X < ^ gives +6 and 
—3 resectively for the vacuum expectation of Ap as the lattice cutoff is removed. 

FIG. 2. Vacuum expectation value of the plaquette, {Ap)^, as a function of X and the lattice 
coupling, ^ with 6 = — 1. In weak coupling, X(. = marks the boundary between two domains of 
convergence. 

FIG. 3. Vacuum expectation value of the plaquette, {Ap)o, as a function of X and the lattice 
coupling, ^ with 6 = 0. In weak coupling, Xc = marks the boundary between two domains of 
convergence. 

FIG. 4. Vacuum expectation value of the plaquette, {Ap)o, as a function of X and the lattice 
coupling, ^ with 6 = 1. In weak coupling, Xc = —\ marks the boundary between two domains of 
convergence. 

FIG. 5. Projectors for the variational basis used in the 0"'"+ glueball calculation. 

FIG. 6. Monte Carlo results for the vacuum expectation values of the seven Wilson loop op- 
erators defining the variational basis. Different curve types correspond to four different values of 
X. 

FIG. 7. SU(3) string tension calculated from Eq.(15) as a function of X (see text) and the 
lattice coupling, ^. All theories corresponding to X > X,, yield the same string tension in the 
continuum limit although the rate of convergence differs. 

FIG. 8. Lowest eigenstate and its ratio with respect to string tension as a function of the lattice 
coupling for X = 0.5, 1.0, 1.5, and 2.0. 
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FIG. 9. Second lowest eigenstate and its ratio with respect to string tension as a function of 
the lattice coupling for X = 0.5, 1.0, 1.5, and 2.0. 

FIG. 10. Third lowest eigenstate and its ratio with respect to string tension as a function of 
the lattice coupling for X = 0.5, 1.0, 1.5, and 2.0. 

FIG. 11. Fourth lowest eigenstate and its ratio with respect to string tension as a function of 
the lattice coupling for X = 0.5, 1.0, 1.5, and 2.0. 
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